Superfluid Density in the s+ Wave state of Clean Iron-based Superconductors 
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Based on a phenomenological model and the Kubo formula, we investigate the superfluid density p s (T) and 
then the penetration depth A(T) of the iron-based superconductors in the coexistence region of the spin-density 
wave and superconductivity, and also in the overdoped region. Our calculations show a dramatic increase of A(0) 
with the decrease of the doping concentration x below x = 0.1. This result is consistent with the experimental 
observations. At low temperatures, p s (T) shows an exponential-law behavior, while at higher temperatures, 
the linear-in-T behavior is dominant before it trends to vanish. It is in qualitative agreement with the direct 
measurement of superfluid density in films of Fe-pnictide superconductor at x = 0.08. The evolution of AA(T) 
can be roughly fitted by a power-law function with the exponent depending on the doping concentration. We 
show that the Uemura relation holds for the iron-based superconductors only at very low doping levels. 

PACS numbers: 74.70.Xa, 74.25.N-, 75.20.-g 



In addition to zero resistance, the Meissner effect is another 
hallmark of superconductivity. The directly measured pene- 
tration depth(/t) in a weak magnetic field provides informa- 
tion of the gap structure, and is a characteristic length scale 
of a bulk superconductor. In general, p s oc I /A 2 . The num- 
ber of electrons in the superconducting phase, p s , character- 
izes the phase rigidity of a superconductor. In conventional 
Bardeen-Cooper-Schrieffer (BCS) superconductors, the pene- 
tration depth exhibits an exponential behavior at low tempera- 
tures, and the power-law behavior in AA{T) = A(T) - A(0) has 
been considered as evidence for unconventional pairing sym- 
metry in the high-temperature superconductors yj . Compare 
to cuprates, the remarkable features of iron pnictides are the 
nature of magnetism and the multiband character. They have 
triggered massive studies since their discovery In this 

letter we focus on its response to a weak external magnetic 
field. 

There are several ways to measure magnetic penetration 
depth 0-0 ■ In the 1111 systems, at low temperatures, some 
experiments J3l found a power-law behavior A(T), while oth- 
ers Jil 01 have found an exponential temperature dependence 
of A(T). The situation in the 122 system is also unclear: 
The superfluid density p s (T) exhibits an exponential behav- 
ior in the cleanest Ba 1 _ x K x Fe2As2 [10], while measurements 
on B a(Fei _ x C o x )? As? have shown a power-law behavior of 
A(T) 111 lWl6ll with the exponent varying from 1.6 to 2.8, and 
a two-gap scenario is suggested for Ba(Fei_ x Co x )2As2 and 
Bai- x Rb x Fe2As2 
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And there are also some theoret- 
ical works II19U22I1. 

In this letter, we carry out systematic calculations of p s (T) 
based on a two-orbital phenomenological model 02311 . Within 
this model, each unit cell accommodates two inequivalent Fe 
ions and results based on this model on various properties of 
Fe-pnictide superconductors l!23M3lll are in reasonable agree- 
ment with experimental measurements. When we normalize 



the energy parameters of the Fe-Fe nearest and next-nearest 
neighbors, the hopping i nteg rals defined below are chosen as 
fi_4 = 1, 0.4, -2.0, 0.04 023J, respectively. In the momentum 
k s pac e, the single-particle Hamiltonian matrix can be written 
as MS 
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with a i = -2?2 cos (k x + k y ) - 2?3 cos (k x - k y ), «2 = 
-2?3 cos (k x - k y ) - 2?2 cos (k x + k y ), 03 = -2?4(cos (k x + k y ) + 
cos (k x - k y )), «4 = -2?i(cos k x + cos ky), where /i is the chem- 
ical potential. Here we have chosen the x axis along the 
link connecting nearest neighbor (NN) Fe ions, and the dis- 
tance between NN Fe is taken as the unit of length. The 
pairing term H Kk = 2»vk( A 

<*M c l n ,ki c av-ki H e.) has only 
next-nearest-neighbor (NNN) intra-orbital pairing, where a 
denotes Fe A or Fe B in the unit cell and v denotes the orbitals. 
It will lead to the s ± -wave pairing symmetry ifiol [Tl , 321. 
The self-consistent conditions are A^k = 2 £ T cos k T A" /+T and 

A ?/ + T = ¥ C lflr. V r C li C lr, V ^ = ^2kCOSk T <C„ v , kT C. v ,_ ki >, 

with t — x ± y and the pairing strength V — 1.2. The inter- 
action term includes the Hund's coupling Jh = 1.3 and the 
on-site Coulomb interaction U, in which we choose U — 3.4 
and U = 4.0 as two different kinds of homogenous systems. 
After taking the mean-field treatment 1124 12511 , //„„ can be ex- 
pressed as 



Hint = U 2^ (*fca>tycr + (U - 3J H ) ^ <n,>r>n/, 

+(U - 2Jh) ( n U&) n 

ivcr • 



(2) 



In the presence of spin-density-wave (SDW) order, in 



2 



the k space can be decoupled into a diagonal term and mag- 
netic term. Define ^ = (4o,kr 4,kr 4o,kr 4i,k T >' fl = 
(^Ir ^k+Qt' ^-ki' 'A-k+Ql)' me Hamiltonian without external 



field in k space can be written as ip'^Hotfk H261 12711 ■ with 
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where / is a 4 x 4 unit matrix, R - -y(C/ + J H )H M , and the 
corresponding //,' k = ff ( k + f(3t/ - 5Jh)I, with « = 2 + x. R 
relates to the magnetic order II261 12711 with 
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in Eq.(4) I is a 2 x 2 unit matrix. Due to SDW order, the 
wave vector k is restricted in the magnetic Brillouin zone 
(BZ). The self-consistent condition is M — \ 2 v (« j 4vt _h Av|) = 
2^ 2 v ,k °" c i V( T-k c ^"7-k+Q, r ab is the distance of Fe B to the ori- 
gin sited by Fe A. N s is the number of unit cells. We take 
N s = 512 to obtain self-consistent parameters and N s = 768 
in the calculation of p s . After diagonalizing 2k ^{flofw - 
2km ^k.mT^Tm by a 16 x 16 canonical transformation matrix 
T, we can obtain all properties of the system without the ex- 
ternal field. 

Our investigation of the superfluid density p* follows the 
linear response approach described by Refs. JTI 33-351- In 
the presence of a slowly varying vector potential A x (r, t) = 
A(q, a))e 1 ^' T ~ 1<ot along the x direction, the hopping term is mod- 
ified by a phase factor, Cj^cj^ — * c\ a Cj a exp i^ f r ' A(r, t) ■ dr. 
Throughout the letter we set ft = c — 1 . By expanding the 
factors to the order of A 2 , we obtained the total Hamiltonian 
H t0 , = H + H' with 



H' = - ^A x (n,t)[eJ^(n) + -e z A x (n,t)K x (ri)]. 



(5) 



J x (rj) is the particle current density along the x axis, K x {ri) is 
the kinetic energy density along the x axis. Their expressions 
are 

Kxin) = - 2^ ki+sxlj+sic^Ci+gyv + H.c), (6) 

vv'crd 

J x( r i) = -i ki+sHt+siclya-Ci+sya- - H.c), (7) 



''<T(j 



only 6 - x,x + y have contributions to the x component and 
Xij+s — 1 in our coordination. The charge current density 
along the x axis is defined as 



J?(rd = 



SH' 



SAJji, t) 



eJ p x (n) + e 2 K x (n)A x (r h t). 



(8) 



The kinetic energy is calculated to zeroth order of A x (rf), cor- 
responding to the diamagnetic part, and that of the paramag- 
netic part J x (ji) is calculated to the first order of AJji). In the 



interaction representation we have 



(J x (n)) 
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(9) 



() represents the expectation value based on the wave func- 
tion of H tot while ()o corresponds to the wave function of 
Ho. In the Matsubara formalism we have the current-current 

rB 

correlation II v . v (q, ia>) = J dTe 10JT T[ xx (q, r), and n o (q, t) = 
-{T T J P X (q,T)7f(-q,0)) = Z mim2 nS m2 (q,r) where T T is the 
time ordering operator, J p (q, r) = e rHa J p x (q)e~ TH °, 7f(q) = 
Si e" lqr '/f in) = 2 mim2 Jm lJn2 (q) is a summation overk. Cal- 
culation of n xl (q, iu>) is in the framework of equations of mo- 
tion of Green's function, 

Tffli/!!2 



dir x r (q,r) 



dj 



= -[^„ m2 (q)^.f(-q)]- 

- <7>"° T [// , ^ lim2 (q)]_ e - H ° T ^(-q, 0)) . 
A lengthy but straightforward algebra leads to 

^k.k+q^+q.k. 



n„(q, iw)=^ 



(f(Ek,nn) ~ /(£k+q,m 2 )) 
ia» + (£k,/m - £k+q,m 2 ) 



(10) 



where / is the Fermi distribution function. Through analytic 
continuation, n xt (q, a>) is obtained. When u> = 0, the deriva- 
tive of / has an important contribution to H xx (q,ico). The 
quantity Y^m? can be expressed as 

Y mim = TT [ f 4(^4(sin + sin £,+,,) + £4 (sin t%_ y + sin fc? +y )) 

+ h{^2 sin fc A -_ v + £2 sin A: x+y + £ 2 sin A:^_ v + ^ 2 sin k® +y ) 
+ ti{^2 sin k x+y + ^2 sin A: v _ y + £ 2 sin A:^ +v + sin A:^_ v ) 
+ fiCfisin^ + fismJ^)], (11) 



with £ 

?2 = 01 

and <? 



k,k+q , k+q,k 

1 = a 13 +a xl - 

k.k+q k,k+q ^ 

3,3 +ff ll,ll'? 4 
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k+q,k k,k+q ^. _ k,k+q k,k+q 

^9,11 11,9 ' '2 — "*" Q; 9 J 9 ' 

k,k+q k+q,k k+q,k k,k+q 

a i,2 + a 2,l tt 9,10 tt 10,9 ' 

TU(k)Wk') + T; +U)li (k)T ;+1 , m2 (k'). The cor- 
responding ^' is connected to ^, by changing a (J into ff, + 4 J+ 4. 
A: A±y denotes A: A + k v and fe^. = k x±y + Q. The superfluid 
weight measures the ratio of the superfluid density to the 
mass DJne 2 = p s /m* = -(J x (ri,t))/e 2 A x (ri), and the Drude 
weight is a measurement of the ratio of density of mobile 
charges to their mass IIU l33l435ll . 

q y ^Q,oj = Q)-(K x )o, (12) 
q y = 0,co^0)-(K x ) . (13) 



Figure Q] shows the variation of D s , D, M and supercon- 
ducting (SC) order A = ± £a(A£ i+x+y + A^ j+X _ y ), as functions 
of x at different temperatures. D does not change much as 
the temperature varies and we plot it clearly in Figs. [Tic) and 
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FIG. 1: (color online) Panels (a), (c), and (d) plot D s (black solid 
line), D (orange dashed line ), A (red dotted line), and M (blue dash- 
dot-dotted line) as functions of x at different temperatures. The right 
scale is for D s and D while the left scale is for A and M. Panel (b) 
plots /l(0) as a function of x. The inset of panel (b) is the phase 
diagram of temperature T and x. 



Old). At zero temperature, we do not show the plot of D be- 
cause in almost all the doping levels D s — D as long as A has 
finite value; Fig.QJa) shows that in the overdoped regime, the 
superconducting gap disappears and D s drops to zero, while 
D is finite just like the plot in panels (c) and (d); hence, in 
the overdoped levels when A = the system corresponds to 
metal. We can see from Fig.QJa) that at T = 0, D s increases 
with the increase of x until it reaches the SDW boundary. In 
the underdoped region x < 0.05, most of the Fermi surfaces 
are gapped by SDW 124 12911 . doping is the major source of 
charge carrier; hence, the superfluid density as well as mobile 
charge density increase linearly with the increase of x. While 
at larger doping 0.5 < x < 0.1, SDW is suppressed, the gapped 
surfaces shrinks significantly, and more intrinsic charge carri- 
ers are released to the system in addition to the doping car- 
riers. This is the reason why the increase of Ds — D with 
doping becomes more dramatic than the linear dependence in 
this region. After SDW disappears, A dominates the behavior 
of D s , and shows a flat behavior in a considerably large doping 
range. In panel (b) we show the variation of A(Q) as a function 
of x for x < 0.3. We define p s (T) = D S (T) = A(TT 2 with 
arbitrary units. Compared to the phase diagram in the inset, 
we find that in the SDW + SC coexisting regime, A(0) shows a 
sharp increase with the decrease of x, which is in good agree- 
ment with experiments 11121 1 1311 . 

An external magnetic field can couple relevant correlation 
functions; hence, p s is a nonlocal quantity, describing the stiff- 
ness of the system. Figure Q2c) and Hid) show that at finite 
T, D s deviates from D, the suppression of D s is stronger than 
that of A. For the U — 4 case, the results (not shown here) are 
very similar to the results presented here. 
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FIG. 2: (color online) Density of states at T = 0.02 for different x. 
All those calculations are for the U = 3.4 case. 
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FIG. 3: (color online) Panels (a), (b), and (c) plot the renormalized 
superfluid density p s (T)/p s (0) and superconducting order parameter 
A(r)/A(0) as functions of the temperature T/T c at different doping 
levels for U = 3.4. T sdw is the transition temperature for SDW. The 
green dotted lines are linear-in-T fitting functions. Panels (a'), (b'), 
and (c') are similar but for V = 4.0. Panel (d),(<i') show the compar- 
ison of our results with experiment data at x = 0.08. Blue solid line 
in the inset of panel (d 1 ) plots p q (T)/ 'p(T c ) as a function of T/T c at 
x = 0.08 and the red dashed line is the aid for the eyes. 



Temperature dependence of superfluid density is a quantity 
reflecting the low-energy residual density of states(DOS) in- 
side the superconducting gap. EquationdTOt indicates that the 
difference between D and D s is related to the derivation of 
/ near the Fermi surface, and can be understood as excita- 
tion of quasiparticles p q . Fig. [2] shows the DOS at T — 0.02. 
For x = 0.05 and 0.1 the gap is considerably larger, hence 
D s is equal or almost equal to D. Although there is a gap at 
x = 0.2[see Fig. |2c)], it is small; therefore, has its 

contribution to D s , and therefore D s deviates from D. 

We choose three typical doping levels, to show the tem- 
perature T/T c dependence of p s {T)/p s {0) and A(T)/A(0) for 
U = 3.4 as well as for U = 4.0. From Fig.|3]we can see that 
the suppression of superfluid density is stronger than that of 
the superconducting order parameter in all cases. At low tem- 
peratures, the curve of p s (T)/p s (Q) is flat, a characteristic of a 
nodeless superconducting gap. 

As T increases, a linear-in-T behavior of superfluid den- 
sity is dominant in all cases. For U = 3.4 cases, linear func- 
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FIG. 4: (color online) Panel (a) plots A/l(T) as a function of T/T c 
at typical selected doping for U = 4, the dashed lines are the corre- 
sponding fitting functions. Panel (b) is the Uemura plot of Fe-base 
superconductor. The x axis is p,,(0) for different doping, the y axis is 
the corresponding T c for the given dopings. 



tions -1.5577r c + 1.52 and - 1. 5777 r c + 1.49 are used to fit 
this kind of behavior for x = 0.1 and x = 0.2, respectively, 
which are shown in Figs. |3b) and |3jc). It is consistent with 
the power-law behavior observed in the experiments 111 lWl6ll . 
Interestingly, they are in good agreement with the direct mea- 
surements of superfluid density in films of Fe-pnictide super- 
conductors in Refll4[ We show our results and the experimen- 
tal data [see Fig. 1(a) in RefuHl together in FigsOJd) (U=3.4 
case) and 3{d') (U=4.0 case), and their consistence is explicit. 
In order to understand the wider linear T dependence of p s {T), 
the inset in Fig.[3jd) plots the renormalized p C] (T)/p(T c ) as a 
function of T/T c at x = 0.08; the red dashed line aids for eyes. 
We can see that the number of excited quasiparticles is expo- 
nentially small at low T with strong superconductivity, but it 
is proportional to linear T within a certain temperature range 
before superconductivity disappears. The easy appearance of 
linear-in-T behavior is closely related to anisotropic S ± super- 
conducting paring, since in-gap states(Andreev states) may be 
induced in this case. The ratio 2At(0) / 'ksl ' c at optimal doping 
is about 4.3 (4.5) for the U = 3.4 (4.0) system. 

Experiments always measure AA(T) = A(T) - A(0), so we 
show the evolution of AA(T) at selected doping concentrations 
for U = 4.0 in Fig.|3a). The results of U — 3.4 are very sim- 
ilar. In the low-temperature range the curve is fiat. At high 
temperature approaching the disappearance of superconduc- 
tivity, there is a jump for the value of AA(T), which we show 
by the colored solid dots. We fit the evolution of AA(T) by 
a power-law behavior. See Fig. EJa); the corresponding fit- 
ting function 4(T /T c ) 3 b (2(T /T c ) 3 ) is for data of x = 0.05 
(x = 0.1,0.2) and it may be the reason why the experiments 
give different exponents for different samples. 



Experiments have shown that the Uemura relation 113 611 
holds 13711 for a 1111 system but does not hold for a 122 
system 0811 . In Fig.|4jb), we plot T c versus p. s (0) based on 
our model. The blue-dashed line (red-dotted line) is for the 
U = 3.4 (U = 4.0) system. It shows that at very low dop- 



ing levels, about x < 0.035(grey point), both the U — 3.4 
and U = 4 systems follow the same empirical linear rela- 
tion(grey line). As T c close to the maximum and p s {0) sat- 
urate at x > 0.08 (0.1) for U = 3.4 (U = 4.0), and the data 
significantly deviate from the linear relation. This is because 
in the very underdoped region the doping is a major source of 
charge carriers and the Uemura relation is valid here. 

Based on a two-orbital phenomenological model, we have 
studied the stiffness of superconductivity in clean iron-based 
superconductors. At zero temperature, we find A(Q) a sharp 
jump as x decreases in the regime of the coexisting SDW + SC 
orders; the variation of A(Q) as a function of doping is in good 
agreement with experiments 01211 . As far as we know this is 
a new theoretical result. At low temperatures, p s (T)/p s (Q) 
is flat, then shows a linear-in-T behavior before the system 
loses its superconductivity. It is in good agreement with ex- 
periments of direct measurement of superfluid density in films 
1 1411 . The evolution of AA(T) roughly follows the power-law 
behavior with different exponents corresponding to different 
doping levels. Only at low doping levels, the empirical Ue- 
mura linear relation holds for the iron-based superconductors. 
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